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Abstract — Evapotranspiration (ET) is the loss of water to 
the atmosphere by combined processes of evaporation 
from soil and plant surfaces and transpiration from 
plants. Since various factors affect ET, including weather, 
crops and soil parameters; numerous equations have 
been developed to quantify standard ET. The equations 
vary in data requirements from very simple, empirically 
based or simplified equations to complex, more 
physically based equations. This study used six methods 
in estimating standard evapotranspiration using data 
from September 201 1 -August 2012 from Climate Station 
at Mas gar (05°10’20” S, 105° 10’ 49 ”E, 50 m dpi) 
Lampung, Indonesia. The six models are: Hargreaves- 
Samani 1985 (H/S), FAO 24 Radiation (24RD), EAO 24 
Blaney-Criddle (24BC), FAO 24 Pan Evaporation 
(24 PAN), Linacre (Lina), and Makkink (Makk). The 
results were analyzed using statistics methods in error 
indicators, which are: Root Mean Square Error(RMSE), 
Mean Absolute Error (MAE), and Logaritmic Root Mean 
Square Error(LOG RMSE), while the closeness among 
the models was analyzed using Index Agreement (IA). 
Direct measurement had been done using lysimeters 
(3x2x1 ) m. The study concluded that Makkink model is the 
suitable simple model that should be chosen in Lampung 
lowland area to calculate ET 0 when climate data is 
limited, besides the recommended EAO 56 Penman 
Monteith. 

Keyword — Evapotranspiration, Standard 

Evapotranspiration, FAO 56 PM, Makkink Model. 

I. INTRODUCTION 

Evaporation is the main component in hydrology cicle, 
therefore, accurate estimation of evaporation rate is 
important for water management and eventually for 
agriculture production. However, it is difficult to 
measure evaporation rate directly since evaporation 
affected by various factors. 

Evaporation is affected by climate factors such as solar 
radiation, air temperature and humidity and wind 
velocity; by crops type and environment and by soil 


condition and management (Temesken, Davidov and 
Frame, 2005). Since those factors are linked to each other 
and change in time and space, it is difficult to develop 
equation for estimating evaporation rate for various crops 
on different condition. Therefore, a scheme called 
reference evapotranspiration was developed. 

Reference ET is defined as “the rate of 
evapotranspiration from an extensive area of 0.08-0.15 
m high, uniform, actively growing, green grass that 
completely shades the soil and is provided with unlimited 
water and nutrients” (Bakhtiari et al., 2011). More 
recently, Allen et al. (1998) elaborated on the concept of 
ET 0 , referring to an ideal 0.12 m high crop with a fixed 
surface resistance of 70 s m' 1 and an albedo of 0.23. 
Several equations has been developed to estimate the 
reference evapotranspiration; some of that were derived 
based on physical processes of the evapotranspiration but 
mostly are empirical based on ststistical relationship 
between evapotranspiration and one or more climate 
variables (Berengena and Gavilan, 2005). Approaching 
methods to estimate evapotranspiration rate was 
developed increasingly in the last 30 years such as based 
on air temperature measurement (Hargraeves and Sumani, 
1985), based on solar radiation Priestly and Taylor 
method (Priestly and Taylor, 1972) and based on 
combination of radiation balance and air moisture 
aerodynamic movement (Penman, 1948). 

Penman method has been improved several times such as 
Penmann method modified by Monteith and known as 
Penman-Monteith equation (Monteith, 1965) approach 
method in FAO 24 version (Doorenbos and Pruitt, 1977), 
FAO 56 modification (Allen, 1998) and recently Matt- 
Shuttle worth approaching method (2009). Those 
estimations 

have been derived and/or calibrated from the direct field 
measurement of ET using various grasses of alfalfas on a 
variety of lysimeter designs, climates, and management 
conditions. 

Some studies showed that Penmann-Monteith model gave 
accurate estimation that FAO and other organizations 
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recommended this model to estimate reference 
evapotranspiration for calculating crops water 
requirement (Itenfisul.et.al., 2003; Berengena dan 
Gavilan, 2005). Even though Penman-Monteith model is 
accepted as accurate estimation, for using in local 
condition, it is necessary to validate the model to whether 
it estimates close enough to the direct measurement on the 
local climate station. Berengena dan Gavilan (2005) 
examined different methods in estimating 
evapotranspiration rate in Southern Spain, an area with 
strong advection. The result showed that Penman method 
modified with local wind function gave the best 
estimation compared with direct measurement in 
lysimeter following with Penman-Monteith FAO 56 
version. Steduto et al. (2003) examine FAO 56 method in 
Southern Italia with Mediteranian semi-arid climate. The 
results showed that FAO 56 method was the best in 
estimating the evapotranspiration; however, tended to be 
over estimated in winter time when the evapotrasnpiration 
was low and under estimated in summer time when the 
evapotraspiration was high. Temesgen et al. (2005) also 
examined FAO 56 in California, USA and the result 
showed good correlations with evapotranspiration rate in 
37 climate stations in this area. 

Direct measurement of evapotranspiration is calculated 
using lysimeters; however, only small number of climate 
stations is equiped with a lysimeter; as an alternative, the 
measurement is done using evaporation pan known as 
Class A pan. A study by Fontenot (2004) showed that 
reference evapotranspiration measurement by Class A pan 
was not fit to the estimation by Penman-Monteith. The 
result from Class A pan should be corrected by a 
coefficient. Using Cuenca and Snyder method; Xing et 
al.(2008) obtained that pan coefficient for Canada maritim 
climate was between 0.78 - 0.94. Generally the result 
from Pan A was lower than the estimation result from 
Penman-Monteith or Priestley-Taylor method. 

Research about comparing different models has been 
done in some countries. Chen et al. (2005) used 7 
estimating models in four provinces of Taiwan and found 
that Makkink and Hargreaves- Samani models were the 
best models in estimating ET 0 when compared to FAO 56 
PM. Chowhury et al. (2010) also found that in India, 
Makkink model had the closest estimation to FAO 56 PM 
with a little underestimated result. 

Xu and Chen (2005) did similar study in Germany with 
comparing 7 models and found that Granger-Gray and 
Makkink models were the best models for the area. In 
North China, Scneider et al. (2007) compared 4 models 
with direct observation and concluded that Hargreaves - 
Samani and Makkink models were the best models in 
estimating ET 0 even better than FAO 56 PM. Jacobs et al. 
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(2004) conducted research on estimating ET 0 in Florida, 
using remote sensing method with data from GOES. The 
results showed that FAO 56 PM is the best model with 
R 2 = 0.92; however, this result was not much different 
with estimated results from Makkink model which gave 
R 2 = 0.90. 

Various methods in estimating evapotranspiration also 
has been applied in Indonesia. Usman (2004) compared 
Thornthwaite, Blaney-Cridle, Samani-Hargreaves, 
Prestley-Taylor, Jansen-Haise, Penman and Penman- 
Monteith methods in five climate stations in West Java; 
the results showed that in general Priestley-Taylor in 
average gave the highest annual evapotranspiration rate, 
while the lowest was obtained by Blaney-Criddle method. 
It also showed that estimation using Penman Monteith 
method in general gave higher rate than Pan A 
measurement. 

Fampung Province (103° 40’ - 105° 50’ E; and between: 
6° 45’ - 3° 45’ S;35.288,35 km 2 ) is located at Southeast 
tip of Sumatra Island, Indonesia. Fampung climate is 
characterized by monsoonal rain distribution and local 
characteristics. Rain season in general is from October to 
March with the peak on January/February and dry season 
is from April to September . Monthly rainfall ranges from 
50 - 200 mm and annual rainfall ranges from 1200 mm 
(lowland area) to 2500 mm (highland area). Fampung 
economic is dominated by agriculture products mainly 
coffee, chocolate, rubber and sugarcane. Fampung is also 
considered as main area for cash crops such as paddy, 
soybean and Maize. Therefore, finding good and reliable 
method in estimating crops water requirement is 
necessary for better agriculture management. 

The objective of this research was to find a closest model 
to the FAO 56 PM model by comparing six different 
models in estimating standard evapotranspiration in 
Fampung area, Indonesia. 

II. METHODS 

This study used six methods in evaluating potential 
evapotranspiration using data from September 2011 to 
Agustus 2012 from Climate Station at Masgar (05°10’20” 
S, 105°10’ 49”E, 50 m dpi) in Fampung, Indonesia. 

The six models are: Hargreaves- Samani 1985 (H/S), FAO 
24 Radiation (24RD), FAO 24 Blaney-Criddle (24BC), 
FAO 24 Pan Evaporation (24PAN), Finacre (Fina), dan 
Makkink (Makk). The results from those models were 
compared to FAO Penman-Monteith (56PM) as the 
standard model. 

To evaluate the relation between models, the results were 
analyzed using statistics methods in error indicators, 
which are: Root Mean Square Error (RMSE), Mean 
Absolute Error (MAE), dan LogaritmicRoot Mean Square 
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Error (LOG RMSE), while the closeness among the 
models was analyzed using Index Agreement (IA). 

Finally, the results were compared to the direct 
measurement using 3 lysimeters (3x2x1 m) planted with 
Sporobulus diander grass. 

2.1 Description of Models 

2.1.1 Hargreaves-Samani 1985 (H/S) (Hargreaves and 
Samani, 1985) 

The equation of this model is: 

ET 0 = 0.0023 (T mean A 17.8)(T max 

Tmin ) 0S R a (l) 

with ET 0 is standard evapotranspiration (mm/day), T mean 
is daily mean temperature (°C), T max is maximum 
temperature(°C), T min is minimum temperature, dan R a is 
radiation on top of the atmosfer (MJ/m 2 /day). 
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—0.0097 [ln{U d + 1 ][ln(RH mln + l] 2 [ln(^+ 1)] 

ET 0 is standard evapotranspiration (mm/day), P is 
percentage of day length, T is daily average temperature 
(°C), RH is minimum relative humidity (%), n/N is ratio 
of possible actual day, Ud is wind speed at 2 m (m/s) 

2.1.4 FAO 24 Pan Evaporation (24PAN) Doorenbos and 
Pruitt (1977) 

The equation of this model is 

ET 0 = Kp X E pan (8) 

= 0.108 - 0.028u 2 + 0.0422 ln(F£T) 

+ 0.1434 In (RH mean ) 

- 0.00063 l[ln(FFT)] 2 In (RH mean ) (9) 


2.1.2 FAO 24 Radiation (24RD) (Doorenbos and Pruitt, 
1977) 

The equation of this model is: 


ET 0 = a + b 


- A 

7 R s 

A Ay s . 


( 2 ) 


ET 0 is standard evapotranspiration (mm/day), A vapor 
pressure curve (kPa/°C),y is psychrometric constant 
(kPa/°C), R s is solar radiation(MJ/m 2 /day), a and b 
conversion factor with a = —0.3 mm/day and b derived 
from the equation: 


ET 0 is standard evapotranspiration (mm/day), K p is pan 
coefficient, E pan is class A Pan evaporation (mm/day), u 2 
is average wind speed (m/s), RH mean relative humidity 
(%), dan FET is distance between pan and green crops 
(m). 


2.1.5 Linacre (LINA) (Linacre, 1977) 


The equation of this model is: 


(SOOT, 


ET 0 = 


Vioo 


5;) + 15 <T - To) 


(80 - T) 


( 10 ) 


b = 1.066 - 0.13 x 10 ~ 2 RH mean + 0.045 U d - 0.20 
x 10 3 RH mean U d 

-0.1315 x 10 ~ 4 RH mean 2 - 0.11 x 10 ~ 2 U d (3) 

RH mean is daily relative humidity (%) andt/ d is average 
wind velocity at 2 m height (m/s) 

2.1.3 FAO 24 Blaney-Criddle (24BC) 

(Jensen et al., 1990) 

The equation for this model is: 

ET 0 = a Abf (4) 

/ = p(0.46T + 8.13) (5) 


T m = T A 0.006/t (11) 

ET 0 is standard evapotranspiration (mm/day) T is mean 
temperature (°C), A is latitude of the climate station (°), 
T m is elevation of climate station (m), dan T d is average 
dew point temperature (°C). T d equation is: 

i 

T d = (j^) 8 (H2 + 0.9 T) + 0.1 T - 112 (12) 

is average dew point temperature (°C), T is mean 
temperature (°C), dan / is average daily relative humidity 
(%)■ 


a = 0.004 RH min - £ - 1.41 (6) 

b = 0.908 - 0.00483fftf min + 0.7949 ^ 
+ 0.768 [ln(U d + l] 2 
-0.0038 RH min ^~ 0.000443 RH min U d 

+ 0.281 [ln(^+l)] (7) 


2.1.6 Makkink (Makk) (Makkink, 1957). 

The equation of this model is: 

et °=°- 6 i A£s-°- u (13) 

ET 0 is average dew point temperature (°C). R s is solar 
radiation (MJ/m 2 /day), A is vapor pressure curve 
(kPa/°C), and y is psychrometric constant (kPa/°C). 
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2.1.7 FAO 56 PM (56PM) (Allen et al., 1998) 
The equation of this model is 


ET 0 = 


0.408A(ft n — G) +y ~~ U 2 
A + y(l + 0.34 U 2 ) 


.(14) 


ET 0 is standard evapotranspiration (mm/day), Rn is netto 
radiation on crops surface (MJ/m 2 /day),G is continuous 
heat flux to soil depth (MJ/m 2 /day),T is daily temperature 
at 2 m (°C),U 2 is wind speed at 2 m (m/s),e s is vapor 
pressure (kPa),e a is actual vapor pressure (kPa), Avapor 
pressure curve (kPa/°C),and y is psychrometric constant 
(kPa/°C). 


In this study the ET 0 estimation from FAO 56 Penman - 
Monteith model as the standard model was calculated 
using CROPWAT. CROPWAT is a ccomputer program 
recommended by FAO based on FAO 56 Penman- 
Monteith model (Allen et al ., 1998) 


2.1.8 Indicators 


The error indicators equation used to evaluate the model 
follows Wilmorth (1982): 


J N 

lY/ET osi - ET omi y (15) 


No 

Model 

Climate data needed by each model 

E pa n 

T 

Rs 

Rn 

RH 

p 

u 2 

Ra 

1 

56PM 


V 

V 

V 

V 

V 

V 

V 

2 

H/S 


V 






V 

3 

24RD 


V 

V 


V 


V 


4 

24 BC 


V 


V 

V 

V 

V 


5 

24PAN 

V 




V 


V 


6 

Makk 


V 

V 






7 

LINA 


V 








Notes: E pa n : Evaporation pan (mm/day); T: average, 
maximum and minimum temperature (°C); R s : solar 
radiation(MJ/m 2 /day); R n : nett radiation (MJ/m 2 /day); 
RH: relative humidity (%); P: day length (%); U 2 : wind 
velocity (m/s); R a : radiation at the top of the atmosphere 
(MJ/m 2 /day) 



MAE 


N 

= ^ ETo -~ ET °^ (16) 

i= 1 


N 

L0G = vZ (log£7osi “ } °z ET °™) 2 ( 17 ) 


Z(ET osi - ET omi y 

n\ET osl '\ + \\ET omi W 


(18) 


ET'osi = ET 0Sl - ET omi (19) 


ETomi ~ ET omi ET omi (20) 


With ET osi is Penman-Monteith standard 
evapotranspiration as the standard model 
-i, dan ET omi is others evapotranspirasi models (i). 

Table. 1: Climate parameters needed by each estimation 

model 



Fig.l: The lysimeters 


III. RESULTS AND DISCUSSIONS 

3.1. Statistical analysis 

The statistical analysis from daily evapotranspiration data 
including the error indicator of each model compare to the 
FAO 56 PM as the standard is presented in Table 2-5. 


Table. 2: RMSE value among the ETo estimating models 


RMSE 


56PM 

Makk 

24BC 

24PAN 

24RD 

H/S 

Lina 

56PM 

0 

0,34 

1,30 

0,75 

0,69 

1,35 

0,88 
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Makk 

0,34 

0 

1,61 

0,48 

0,49 

1,52 

1,12 

24BC 

1,30 

1,61 

0 

1,99 

1,92 

1,12 

0,79 

24PAN 

0,75 

0,48 

1,99 

0 

0,33 

1,93 

1,54 

24RD 

0,69 

0,49 

1,92 

0,33 

0 

1,98 

1,54 

H/S 

1,35 

1,52 

1,12 

1,93 

1,98 

0 

0,59 

Lina 

0,88 

1,12 

0,79 

1,54 

1,54 

0,59 

0 


Comparison of error indicator (RMSE) between models using 
monthly data was presented in Table 2. Based on the 
comparison among the six models, the error indicator RMSE 
ranged from 0.33-1.99 which means that ET 0 difference among 
the models was 0.32 mm to 1.99 mm/day. This is not a small 
number since 1 mm/day ET in 1 ha area is equivalent with 
water loss of 10,000 liter/day or 3.6 million liter/year. 

Using Lampung climate data, the lowest RMSE was found 
between FAO 24 Radiation and FAO 24 Pan Evaporation 
while the highest RMSE was found between model FAO 24 
Pan Evaporation and FAO 24 Blaney-Criddle. 

Tabel.3: MAE value of estimated monthly evapo- transpiration 
data among the models 


MAE 


56P 

M 

Mak 

k 

24B 

C 

24PA 

N 

24R 

D 

H/S 

Lin 

a 

56PM 

0 

0,28 

1,06 

0,62 

0,67 

1,28 

0,8 

6 

Makk 

0,28 

0 

1,28 

0,40 

0,45 

1,50 

1,0 

8 

24BC 

1,06 

1,28 

0 

1,69 

1,74 

0,95 

0,6 

8 

24PA 

N 

0,62 

0,40 

1,69 

0 

0,25 

1,90 

1,4 

8 


24RD 

0,67 

0,45 

1,74 

0,25 

0 

1,95 

1,5 

3 

H/S 

1,28 

1,50 

0,95 

1,90 

1,95 

0 

0,4 

7 

Lina 

0,86 

1,08 

0,68 

1,48 

1,53 

0,47 

0 


The second error indicator (MAE) is presented in Table 3. 
Similar results with RMSE were found in error indicators 
for both MAE and log RMSE (Table 4). Makkink model 
was the model which is closest to FAO 56 Penman 
Monteith. 


Table. 4: LOG RMSE among the estimating models ofETo 


LOG RMSE 


56PM 

Makk 

24BC 

24PAN 

24RD 

56PM 

0 

0,04 

0,13 

0,11 

0,10 

Makk 

0,04 

0 

0,17 

0,08 

0,07 

24BC 

0,13 

0,17 

0 

0,23 

0,22 

24PAN 

0,11 

0,08 

0,23 

0 

0,06 

24RD 

0,10 

0,07 

0,22 

0,06 

0 

H/S 

0,14 

0,16 

0,11 

0,23 

0,23 

Lina 

0,10 

0,13 

0,08 

0,20 

0,19 


MAE between FAO 56 Penman-Monteith and other 
models ranges from 0.28 mm/day (Makkink) up to 1.28 
mm/day (Hargreaves-Samani 1985) and LOG RMSE 
ranges from 0.04 mm/day (Makkink) to 0.14mm/day 
(Hargreaves-Samani). 

Table 5 showed the results of Index of Agreement (IA). 
Consistently, Makkink model gave the best results with 
IA 0.78 followed by Linarch (0.42) and FAO Pan 
Evaporation (0.42) 


Tabel.5: Index of Agreement among the models 


I.A 


PM 

MK 

BC 

Pan 

24 RD 

HS 

Ln 

PM 

1 

0,78 

0,09 

0,42 

0,55 

0,26 

0,42 

MK 

0,78 

1 

-0,35 

0,80 

0,81 

0,10 

0,10 

BC 

0,09 

-0,35 

1 

-0,40 

-0,22 

0,79 

0,85 

Pan 

0,42 

0,80 

-0,40 

1 

0,95 

-0,03 

-0,08 

24 RD 

0,55 

0,81 

-0,22 

0,95 

1 

-0,03 

-0,01 

HS 

0,26 

0,10 

0,79 

-0,03 

-0,03 

1 

0,93 

Ln 

0,42 

0,10 

0,85 

-0,08 

-0,01 

0,93 

1 


The resume of statistical analysis from daily evapo transpiration data including the error indicator of each model compare to 
the FAO 56 PM as the standard was presented in Table 6. 
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Table. 6: Statistical of daily evapotranspiration data from each model 


Evapotranspiration Model 

Average 

STD 

RMSE 

MAE 

LOG RMSE 

I.A 

24 

BC 

HS 

Mk 

Lin 

Pan 

FAO 

2,851 

4,607 

4,821 

3,306 

4,387 

2,925 

3,533 

0,485 

1,419 

0,561 

0,370 

0,569 

1,009 

0,774 

0,884 

1,392 

1,481 

0,683 

1,065 

1,283 


0,749 

1,107 

1,294 

0,571 

0,890 

1,037 


0,117 

0,130 

0,170 

0,088 

0,131 

0,239 


0,439 

0,565 

-0,403 

0,595 

0,263 

0,363 



Table 6 showed that among the models, Makink model 
consistantly had the smallest RMSE, MAE and Log 
RMSE compared to PM model and had the highest 
agreement. On the other hand HS model had the biggest 
RMSE, MAE and log RMSE with the lowest agreement. 
Therefore, for Lampung, estimation ET model with the 
closest estimation to FAO 56 Penman-Monteith was 
Makkink model. It can be concluded that Makkink model 
was the suitable simple model that should be chosen in 
Lampung to calculate ET 0 besides the recommended one, 
FAO 56 Penman Monteith, especially when the climate 
data is limited. 

So far the estimating model that broadly used is FAO 24 
PAN which is based on observation on Pan A evaporation 
pan. In this study this model did not give a good 
estimation compared to the FAO 56 PM model (RMSE 
0.75; MAE 0.62; Log RMSE 0.11 and IA 0.42). In 
comparing 24 PAN model to 56 PM, using 3 years data in 
2 stations in Lampung, Manik et al. (2012) found that the 
coefficient correlation between those two models are low 
(r=0.3 for Branti Station and 0.5 for Masgar station). 

Daily modeled ET 0 results from each model in 1 year is 
presented in Figure 2. Most of the models had similar 
trends with FAO 56 PM but with different estimation. 
Some models underestimated FAO 56 PM (Makkink, 
FAO 24 Radiation and FAO 24 Pan Evaporation) while 
some overestimated ( Blanney Criddle, Hargreaves - 
Samani 1985 and Linarch). Makkink model had good 
estimation to FAO 56 PM in wet months October-March, 
and slightly underestimated in dry months March - 
October. 


* 

s 

1 r 

E 5 
o 4 

| 

\ 

0 


Fig. 2: Daily estimation of evapotranspiration from 
different models 


——24 EC — HS f/k — Ut —FAO 



(Uyt 


In general, Makkink model is a function of R s (solar 
radiation MJ/m 2 /day), A(slope of vapour pressure), and 
y(psychrometric constant). Makkink is a simple model 
since y is 66,1 (kPa/°C), while R s and A could be 
calculated using following equations (Allen, et al., 1998): 


Rs — Krss/ O' max Tmin)Ra ( 21 ) 


4098 [0,6108 exp 

A _ Y_f_ r \T +237,3/ ] 

(T + 237, 3) 2 1 j 

K Rs is a coefficient, 0.16 for interior land area and 0.19 for 
coastal area, R a is top solar radiation (MJ/m 2 /day). 
Basically, this model can be calculated using only 
maximum and minimum temperature (°C) which is more 
available in most research stations. 

Irmak, Allen and Whitty (2003) conducted a research 
using daily measured weather data for a 23 -year (1978- 
2000) in North-Central Florida to examine twenty one 
ETo methods (excluding the FAO 56-PM; 10 combination 
methods, 4 radiation methods, 5 temperature methods and 
2 pan evaporation methods) and the results showed that 
the performance of all radiation methods including 
Makkink was poorer than that of all combination methods 
except the Stephens-Stewart method, which performed 
better than the original PM combination method. 
Makkink methods had a similar standard error of estimate 
for all months with the Stephens-Stewart method, but the 
method significantly underestimated ETo throughout the 
year. 

The tendency to underestimate high evaporative demand 
in Makkink model has been found in most of previous 
studies e.g. in Korea (Chen, et al., 2005), in Germany (Xu 
and Chen, 2005) in Jordan (Mohawesh, 2011) in Iran 
(Bakhtiari et al., 2011); those results might be related to 
the ignorance of the significant influence of wind speed 
on ETo. Regardless of that, Makkink model was 
considered as a good option model and the closest to 
Penman Monteith method in India (Haidar, Kumar and 
Sehgal, 2005) and in Hungary (Racz, et al., 2013). 
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3.2. Comparing Makkink model with direct measurement 
Direct evapotranspiration measurement was done using 
lysimeters. Measuring water input and out put together 
with soil water content on the lysimeters was a challenge. 
During this research the measurements were repeated 
several times, however due to technical problems they 
were not always done in the same day. Therefore, the 
results (Table 7 and 8) were considered as an average 
number. 
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In average from Table 7 and 8, evapotranspiration rate 
during the research was 3.8+1.11 =4.91 mm/day, higher 
than estimated Makkink (3.306 mm) and FAO (3.533 
mm) (Table 6). In general, the average and accumulation 
of observed Pan A gave slightly higher 
evapotranspirasion compared to estimated Makkink. 
However, observed Pan A had much wider standard 
deviation; Makkink model gave more flat result in daily 
estimation while Pan A was more fluctuative. 


Table. 7: Water input and output to the lysimeters 


Time of Observation 

(n) 

Rainfal (R) 
(mm) 

Irrigation (I) 
(mm) 

Percolation(P) 

(mm) 

R+I-P 

(mm) 

1 

0.00 

0.00 

0.00 

0.00 

2 

0.00 

0.00 

0.00 

0.00 

3 

0.00 

0.00 

0.00 

0.00 

4 

0.00 

0.00 

0.00 

0.00 

5 

3.30 

0.00 

0.00 

0.00 

6 

0.00 

0.00 

0.00 

0.00 

7 

43.00 

0.00 

0.00 

0.00 

8 

9.40 

0.00 

0.00 

0.00 

9 

26.50 

0.00 

0.00 

0.00 

10 

1.00 

0.00 

0.00 

0.00 

11 

12.50 

0.00 

2.70 

9.80 

12 

0.00 

0.00 

4.17 

-4.17 

13 

0.00 

0.00 

1.62 

-1.62 

14 

0.00 

0.00 

0.76 

-0.76 

15 

0.00 

10.00 

0.38 

9.62 

16 

0.00 

15.00 

0.24 

14.76 

17 

2.00 

0.00 

1.21 

0.79 

18 

9.00 

0.00 

1.19 

7.82 

19 

0.00 

0.00 

2.07 

-2.07 

20 

0.00 

0.00 

0.66 

-0.66 

21 

0.00 

0.00 

0.26 

-0.26 

22 

5.00 

0.00 

0.00 

0.00 

23 

2.70 

0.00 

0.18 

2.52 

24 

5.00 

0.00 

0.03 

4.97 

25 

0.50 

0.00 

0.02 

0.48 

26 

0.00 

15.00 

0.00 

15.00 

27 

0.00 

0.00 

0.00 

0.00 

28 

0.00 

10.00 

0.00 

10.00 

29 

3.00 

0.00 

2.30 

0.70 

30 

3.50 

0.00 

1.75 

1.75 

31 

0.00 

0.00 

0.00 

0.00 

32 

5.50 

0.00 

0.00 

5.50 

33 

0.00 

0.00 

0.23 

-0.23 
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Time of Observation 

(n) 

Rainfal (R) 
(mm) 

Irrigation (I) 
(mm) 

Percolation(P) 

(mm) 

R+I-P 

(mm) 

34 

3.40 

0.00 

0.00 

3.40 

35 

0.50 

0.00 

0.17 

0.33 

36 

3.20 

0.00 

0.02 

3.18 

37 

0.00 

11.67 

0.02 

11.65 

38 

0.00 

53.33 

0.00 

53.33 

39 

0.00 

0.00 

3.00 

-3.00 

40 

2.60 

0.00 

1.53 

1.08 

41 

0.50 

0.00 

0.73 

-0.23 

42 

3.20 

0.00 

0.00 

0.00 

43 

0.00 

0.00 

0.44 

-0.44 

44 

0.00 

0.00 

0.00 

0.00 

45 

0.00 

0.00 

0.08 

-0.08 

46 

0.00 

0.00 

0.11 

-0.11 

47 

24.00 

0.00 

0.00 

24.00 

48 

0.00 

0.00 

0.71 

-0.71 

49 

0.00 

0.00 

0.00 

0.00 

50 

0.60 

0.00 

0.00 

0.60 

51 

0.00 

0.00 

0.00 

0.00 

52 

0.00 

0.00 

0.00 

0.00 

53 

1.00 

0.00 

0.00 

1.00 

54 

0.00 

0.00 

0.00 

0.00 

55 

29.50 

0.00 

0.00 

29.50 

56 

44.90 

0.00 

0.00 

44.9 

57 

0.00 

0.00 

0.00 

0.00 

58 

0.00 

0.00 

0.00 

0.00 

59 

22.00 

0.00 

0.00 

22.00 

60 

8.00 

0.00 

0.00 

8.00 

61 

0.40 

0.00 

0.00 

0.40 

Average 

3.80 


Table. 8: Soil water content of the lysimeters 


Number of 

Observatio 

n 

Soil 

Tension 

Volumetri 

c 

Soil water 

Daily 

AS 

Cumulativ 

e 

AS 

Number of 

Observatio 

n 

Soil 

Tension 

Volumetri 

c 

Soil water 

Daily 

AS 

Cumulativ 

e 

AS 



content 





content 



(n) 

(kQ) 

(%) 

(mm) 

(mm) 

(n) 

(kQ) 

(%) 

(mm) 

(mm) 

1 

122 

43.96 



13 

171 

39.20 

-0.41 

-11.90 

2 

113 

45.04 

2.70 

2.70 

14 

168.33 

39.42 

0.55 

-11.35 

3 

131 

42.96 

-5.21 

-2.51 

15 

179 

38.56 

-2.17 

-13.51 

4 

143.33 

41.69 

-3.17 

-5.68 

16 

176.67 

38.74 

0.46 

-13.05 

5 

146.67 

41.37 

-0.81 

-6.49 

17 

188.33 

37.84 

-2.25 

-15.30 

6 

157.33 

40.38 

-2.47 

-8.96 

18 

203.33 

36.76 

-2.70 

-18.01 

7 

160 

40.14 

-0.59 

-9.56 

19 

236.67 

34.62 

-5.35 

-23.36 

8 

154 

40.68 

1.35 

-8.21 

20 

280 

32.25 

-5.93 

-29.28 

9 

153.67 

40.71 

0.08 

-8.14 

21 

246.67 

34.04 

4.47 

-24.82 

10 

152.67 

40.80 

0.23 

-7.91 

22 

249 

33.90 

-0.33 

-25.15 

11 

148.33 

41.21 

1.02 

-6.89 

23 

251.67 

33.75 

-0.38 

-25.52 

12 

169 

39.37 

-4.60 

-11.49 




average 

-1.11 
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B. Using Pan A evaporation 


Table. 8: Estimation of ETo comparing Pan A with The Makink 



April 

May 

June 


Pan A 

Makink 

Pan A 

Makink 

Pan A 

Makink 

Average Sum Standard 







Deviation 

3.96 

3.39 

3.40 

3.45 

4.19 

3.40 


118.72 

101.80 

105.50 

106.99 

121.38 

98.66 


2.15 

0.25 

1.61 

0.36 

2.01 

0.30 


April 2013 
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Fig. 3: Daily evapotranspiration from observed Pan A compared to estimated Makkink. 


IV. CONCLUSION 

This study concluded that Makkink model is a simple 
model that can be chosen in Lampung as an alternative to 
estimate standard evapotranspiration in an area with 
limited climate data needed to apply FAO 56 PM, with a 
note that Makkink tended to be underestimated during dry 
months. Estimation of evapotranspiration using models 
are sufficient for averaged and accumulated result from 
some period of time, not for daily or single measuremnt. 
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